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Abstract 

Phenomenological models of heavy flavour decays differ significantly in their 
predictions of global features of B c decays, like the B c lifetime or the relative 
weight of c — > s and b — > c transitions. The l/mq expansion which is directly 
based on QCD allows predictions on the pattern to be expected, namely t(B c ) 
to lie well below 1 psec with c — * s dominating over b — > c and a reduced 
semileptonic branching ratio. Due to interference effects one also predicts a 
lower charm content in the final states of B c decays than naively anticipated. 
The numerical aspect of the predictions, however, has to be viewed with consid- 
erable caution since one cannot expect the l/m c expansion to converge readily 
for AC = 1 transitions. 



1 Introduction 



B c mesons consisting of two heavy quarks - B c — (be) - are not easily produced. 
On the other hand it is highly desirable to obtain large samples of them. For their 
study would deepen our quantitative understanding of the inner workings of QCD in 
a significant way: one expects a rich spectroscopy for the (be) boundstates probing 
the inter-quark potential at distances intermediate to those determining quarkonia 
spectroscopy in the charm and beauty systems [|TJ; the Isgur-Wise function for the 
striking channel B c — > Ivip can be calculated; the weak B c decays reflect a multi- 
faceted interplay of various dynamical mechanisms. It is this last aspect I will analyze 
in this note: in Sect. 2 I will review previous phenomenological descriptions of B c 
decays; in Sect. 3 I introduce a genuine QCD treatment based on the heavy quark 
expansion and apply it to inclusive B c transitions; in Sect. 4 I discuss the final states 
in B c decays before presenting my conclusions in Sect. 5. 
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2 A First Phenomenological Look at B c Decays 



There are three classes of transitions contributing to B c decays with roughly compa- 
rable strength; they appear to be easily distinguishable on the diagrammatic level. 
The first two are the decay of the b quark and that of the c antiquark. Since b — > c 
and c — > s transitions do not interfere with each other in any appreciable way, one 
can cleanly separate their widths. The only subtlety here is that b — > ccs decays lead 
to two c antiquarks in B c decays; there arises then an interference between different 
decay amplitudes that is usually referred to as PI. The third class of transitions is 
produced by Weak Annihilation (WA) of b with c. To lowest order in the strong in- 
teractions the WA amplitude suffers helicity and wavefunction suppression (the latter 
reflecting the practically zero range of the weak interactions). Yet for B c — ► sc they 
are represented by m c /m b and f(B c )/mb and thus relatively mild (f(B c ) ~ 450 — 700 
MeV |1[]); furthermore these reductions are partially offset by the factor 16ir 2 reflecting 
the enhancement of two-body phase space - relevant for WA - over three-body phase 
space appropriate for the spectator decay. As explained later, interference between 
b decay and WA can arise; yet this is usually ignored in phenomenological analyses. 
One thus writes down for the nonleptonic width: 



t nl (b c ) = r;^(5 c ) + Y s b z% d m + r££,(B c ) - \at^ c5s (b c )\ + f wa (b c ) (i) 



The expression simplifies for the semileptonic width since no interference occurs there 



Very naively one might equate r g _» s (.B c ) with T(D°) and T b ^ c (B c ) with r(5<j) and 
thus 



and B c decays would be dominated by c — > s over b — > c in the ratio of roughly 4:1. It 
is quite natural, though, to suspect that eq.(3) represents a gross oversimplification. 
Two specific alternatives have been suggested: 

(i) The phase space in B c = (be) — > bsdu ~ B s (du) is more limited than in D = 
(qc) — > qsdu ~ K{du). This could - due to the high sensitivity of the c decay width 
to the available phase space - reduce T 5 ^s(B c ) significantly relative to T(D). No such 
reduction is expected for r b ^ c (B c ). Therefore the relative weight of b — > c transitions 

Only B c — * t + v leading to a purely leptonic final state possesses an appreciable rate. 




(2) 




(3) 



t(B c ) ~ 3 • 10~ 13 sec 



(4) 
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in B c decays gets enhanced. Using a simple recipe for estimating the phase space 
dependance of the quark decay width the authors of ref.@ estimate 

t(B c ) ~ 5 ■ KT 13 sec (5) 

with - and that is the major difference to the naive guestimate given above - c — > s 
transitions now holding only a slight edge over b — > c decays. 

(ii) It has been advocated by Eichten and Quigg that in the expression for the decay 
width - Tq oc Gpm% - one should use a quark mass reduced by the binding energy 
inside the hadron. For the ordinary mesons B^-, B u and B s this can be seen effectively 
as a redefinition of the quark mass. Yet for the more tightly bound system B c there 
arises an observable difference: the binding energy fi B E being the same for the charm 
and for the beauty mass represents a larger fraction of the charm mass than of the 
beauty mass: (m c - li BE )/m c = 1 - n BE /m c < (m b - fJL B E)/m b = 1 - ^ BE /m b . 
Inside the B c meson the c — > s rate will therefore be more reduced than the h — > c 
rate. Furthermore the impact of this mass shift on the width is greatly enhanced: 
ATq/Tq ~ 5 • [ibe/itiq due to Tq oc rrigl For /jlbe = 500 MeV one finds T 5 ^ s and 
T^c reduced by a factor of 6 and 1.7, respectively! This leads to the guestimate 

t(B c ) ~ 1.3 psec , (6) 

i.e., a considerably longer B c lifetime; the b —>■ c transitions now occur somewhat 
more frequently than the c — > s ones. 

The two questions raised above - (i) whether the B c lifetimes is short, i.e. well 
below 1 psec, or 'long', i.e. roughly 1 psec or longer, and (ii) whether B c decays are 
driven mainly by b — > c or by c —>■ s transitions - are highly important and deserve 
study with the best available theoretical tool, the I/tuq expansion. 

3 Treating B c Decays through a 1/rag Expansion 
3.1 General Methodology 

In analogy to the treatment of e + e~ — > hadrons one describes the transition rate into 
an inclusive final state / through the imaginary part of a forward scattering operator 
evaluated to second order in the weak interactions 0, [3J]: 

f(Q^ f -+Q)=ilm J #x{C w {x)tf w (G)} T (7) 

where {.}t denotes the time ordered product and Cw the relevant effective weak 
Lagrangian expressed on the parton level. If the energy release in the decay is suf- 
ficiently large one can express the non-local operator product in eq.(7) as an infinite 
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sum of local operators Oj of increasing dimension with coefficients containing higher 
and higher powers of I/tuq. The width for Hq — > f is obtained by taking the ex- 
pectation value of T between the state Hq. For semileptonic and nonleptonic decays 
treated through order l/m^ one arrives at the following generic expression 



T(H Q /) = ^^-\KM\ 2 
v Q J J 1927T 3 1 1 



f/M \nn\M \m f (H Q \Qio-GQ\H Q ) 

C3{Hq\QQ\H Q } + c J 5 — 2 



+ ^ 4 }^\m q )f,Q)\H Q ) + % 



(8) 



where the dimensionless coefficients c{ depend on the parton level characteristics 
of / (such as the ratios of the final-state quark masses to wq); KM denotes the 
appropriate combination of KM parameters, and a • G — o^G^ with G^ u being the 
gluonic field strength tensor. The last term in eq.(8) implies also the summation over 
the four-fermion operators with different light flavours q. It is through the quantities 
(HQ\Oi\Ho) that the dependence on the decaying hadron Hq, and on non-perturbative 
forces in general, enters; they reflect the fact that the weak decay of the heavy quark 
Q does not proceed in empty space, but within a cloud of light degrees of freedom - 
(anti) quarks and gluons - with which Q and its decay products can interact strongly. 
These are matrix elements for on-shell hadrons Hq; T(Hq — > /) is thus expanded 
into a power series in Hhad/ m Q < 1- For rriQ — > oo the contribution from the lowest 
dimensional operator obviously dominates; here it is the dimension-three operator 
QQ. Since (Hq\QQ\Hq) = 1 + C(l/mg) holds, one reads off from eq.(8) that the 
leading contribution to the total decay width is universal for all hadrons of a given 
heavy-flavour quantum number; i.e., for itlq — > oo one has derived - from QCD proper 
- the spectator picture. Contributions from what is referred to as WA and PI in the 
original phenomenological descriptions are systematically and consistently included 
through the dimension-six four-fermion operators in eq.(8). 

Yet the I/tuq expansion goes well beyond reproducing familiar results. It shows 
the leading nonperturbative corrections to integrated inclusive rates to arise in order 
I/tuq controlled by the expectation values of dimension- five operators [|J. These 
terms had been overlooked before. What is crucial for our subsequent analysis is 
the absence of contributions of order 1/rriQ. This is due to the fact that there is no 
relevant dimension-four operator that cannot be removed by applying the equation 
of motion [J5], |3j; it can also be understood as due to a subtle intervention of the 
local colour gauge symmetry. A phenomenological ansatz on the other hand where 
the quark mass appearing in the decay width is reduced by a 'binding energy' leads 
to large corrections of order l/m^; see the discussion above eq.(6). This is in clear 
conflict with what holds in QCD! 

Using the equation of motion one can obtain a I/ttlq expansion also for the leading 
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operator in eq.(8), QQ. Its expectation value can then be expressed as follows 0: 

(H Q \QQ\H Q ) = 1 - + M£*a + 0(1/^8) ( 9 ) 



where ((pq) 2 )h q = {HQ\Q{iD) 2 Q\Ho) denotes the average kinetic energy of the quark 
Q moving inside the hadron Hq and {^ g )h q = {Hq\Q^u ■ GQ\Hq). 

Eqs.(8,9) show that the nonperturbative contributions to the width through order 
1/uiq can be expressed through three expectation values: (^g)h q , ((pq) 2 )h q and 
(Hq(p)\Ql1^<1l)((1l1uQl)\Hq(p)). The size of the mesonic matrix element of the 
chromomagnetic operator is obtained from the hyper-fine splitting: 

( M 2 g ) Pq ~^(M^-M| q ), (10) 

where Pq and Vq denote the pseudoscalar and vector mesons, respectively. For the 
average kinetic energy we have the model-independant bound [[J 

{(Pq?)h q > (» 2 g ) Hq (11) 

and we know it cannot be much larger than that. The I/ttlq corrections thus largely 
cancel in the expectation value of the operator QQ between meson states: 

(P Q \QQ\P Q )^l + 0(l/m 3 Q ) 

The expectation values for the four-quark operators taken between meson states can 
be expressed in terms of a single quantity, namely the decay constant: 

(H Q {p)\Q Ll ,q L ){q Llv Q L )\H Q {p)) ~ \f HqPii p v (12) 
where factorization has been assumed. 



3.2 B c Decays through Order l/rriq 

Since b — » c and c — > s decays do no interfere with each other in any practical way, 
one can cleanly separate their widths: 

r(fle) = r£?W + r£?(2J e ) + o(i/r^j (is) 

These widths are denoted as Y decay rather than Y spect for a reason: they describe the 
quark decays as proceeding in an environment shaped by the other components of the 

2 I use the relativistic normalization for the states. 
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decaying hadron Hq as expressed by the expectation values {{j>q) 2 )h q and {[1g)h q ', 
thus they go beyond the simple spectator picture. 

The decay widths include 1 / itlq corrections which consist of the semileptonic and 
nonleptonic components: 



^decay 



b^cl 



-pdecay / R 
1 b^cqq'\ n c 



(B c \bb\B c ) 



J o (x c ,0,0) + 



hit. 



x—-2)I (x c ,0,0) 



dx 



(14a) 



N c ■ (B c \bb\B c ) 



A [J2lo(x c ) + 



(/4> 



mi 



dx 



-2)£j (s c )]- 



8A 



mi 



■ [I 2 (x c ,0,0,)+I 2 (x c ,x c ,0)} 



(146) 



(14c) 



where the following notations have been used: Jo an d I 2 are phase-space factors: 

I (x, 0, 0) = (1 - x 2 )(l - 8x + x 2 ) - 12x 2 log a; (15a) 
/ 2 (x,0,0) = (1 -x) 3 , x c = {m c /m b f (156) 



I (x, x, 0) = v(l - Ux - 2x 2 - 12a; 3 ) + 24x 2 (l - x 2 ) loj 



1 + v 



X r, x . n. , 1 + V 



7 2 (a;, x, 0) = v(l H h 3a; 2 ) - 3a; (1 - 2a; 2 ) log 

2 1 — u 



Vl -4x (15c) 
(15d) 



with 10,2(^,^,0) describing the 6 — > ccs transition, and J2Io(%) = Io(%, 0,0) + 



l (x,x,0); A = 7] J, A 2 



c 2 )/6, where 77 = (c 2 . + 2c+)/3, and J represents the 



effect of the subleading logarithms 0. With x c ~ 0.08 one obtains 



/o(a;,0,0)| a;= o.o8-0.56 , l 2 (x,0 
I {x,x, 0)1^=0.08 - 0.24 , I 2 (x,x, 



,0)1^=0.08 ^ 0.78 for b -> cud 
U=o.os - 0.32 /or 6 > ccs . 



Since these functions are normalized to unity for a; = 0, one notes that the final-state 
quark masses reduce the available phase space quite considerably in this reaction. 
The expressions are simpler for c — > s: 



rSXW = IT • (^IcciSe) 



(<■) 



l (x s ,0,0) + 



(Pg)b b , d 



m 2 „ dx 



x— - 2)I (x s , 0, 0) 



(16a) 



r ™z(£ C ) = r o Cj " Nc ■ (B c \cc\B c ) 



A [/ (a; s ,0,0) + 



%(xA-2)J (x Sl O,0)]. 
m; aa; 
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</4> 



• 7 2 (a; s ,0,0, 



(166) 



\92tx 6 m l c 

and the radiative corrections lumped into Ao and A2 are given by the appropriate 
values for c + and c_. With x s ~ 0.012 one finds: 

/ (z,0,0)U=o.oi2^0.91 , / 2 (a;,0,0)| x= o.oi2^0.96, 
i.e. there is much less phase space suppression than for b — > c transitions. 

The transition operators driving £? c decays are the same that generate B and 
D decays. However their expectation values are evaluated for the B c state, rather 
than the B and D state reflecting that the b — > c and c — > s transitions proceed 
in a different environment. The expectation value of the chromomagnetic operator 
is again given by hyperfine splitting between the masses of B* and B c . Those have 
not been measured yet; on the other hand the theoretical predictions should be quite 
reliable for those. With M(B*) ~ 6.33 GeV and M(B C ) ~ 6.25 GeV one obtains 

(B c \b l -a ■ Gb\B c ) ~ (B c \c l -a ■ Gc\B c ) ~ 0.75 (GeV) 2 , (17) 

which is twice the value as for mesons with light antiquarks: 

(B\b-a ■ Gb\B) -0.37 (GeV) 2 , (D\c-a ■ Gc\D) ~ 0.41 (GeV) 2 . (18) 

(B c \Ha ■ Gb\B c ) / 

2 — — ^ 0.033 (19a) 



Thus 



mi 

(B c \c\o ■ Gc\B c ) 



ml 



0.38 ; (196) 



i.e., this correction becomes quite large in the c — > s transition. 

Putting everything together one finds for the B c width through order l/m 2 c : 

T(B C ) ~ 0.95 • r(5 d ) + 0.75 • T(D°) + 0(\jm\ c ) ~ (4.1 • 10" 13 sec)" 1 (20) 

%^~0,6; (21) 

i.e., a short lifetime with c — > s transitions dominating all S c decays! One also obtains 
a rather low semileptonic branching ratio 

BR SL (B C ) ~ 6 % (22) 

with half oi the semileptonic -B c decays being generated by 6 — > c/z/. However these 
numbers have to be taken with quite a grain of salt. For the nonperturbative correc- 
tions in the c — > s component of the Z? c width are very large, as indicated by eq.(19b). 
Thus one can hope only for a semi-quantitative treatment of that component. 



3.3 Order l/rriq Contributions 

In order 1/m.g the explicitely flavour dependant terms appear that had been antic- 
ipated in previous phenomenological studies. Due to the large value predicted for 
f(B c ) they are quite sizable: as indicated in eq.(l) PI reduces the rate for b — > ccs 
to proceed inside B c mesons by ~ 20 - 40% and WA contributes in an only mildly 
suppressed manner. In addition a more subtle effect arises that had not been incor- 
porated into eq.(l): T W a(B c ) and Tl^,f(B c ) can no longer be separated in a strict 
manner. For those two classes of reactions - for b —>■ cud as well as for b —>■ ccs 
modes - can interfere with each other once gluon emission generates a cc pair in WA: 
B c — > bcg*g* — > sc(cc) g * g *. While this observation has no relevant impact on the 
predicted overall lifetime ||, it becomes very important in the analysis of the final 
states to be given in the next section. 

As far as the total lifetime is concerned, the most relevant effect is produced by 
WA which to lowest order in the strong interactions leads to a moderate reduction in 
lifetime: 

t(B c ) ~ 4 • 10~ 13 sec . (23) 

4 On the Pattern in the Final States 

As stated already in eq.(l) PI reduces T b ^ c£s (B c ) considerably. The interference 
between WA and b decays sketched above also reduces the charm content in the final 
state of AB = 1 B c decays in general. The argument goes a follows: To describe the 
impact of WA on the decay rate beyond the lowest order in the strong interactions one 
has to include the emission of 'off-shell' as well as 'on-shell' gluons with the former 
hadronizing into a cc pair: 

be — > duj sc + g*g* — > duj sc + cc 

For this reaction can interfere with the lowest order decay process. As shown in ref.||, 
this interference is destructive and it actually will reduce the rate for B c — > sccc and 
quite possibly also for B c — > dusc by roughly 5-10%. This decrease in the decay rate 
is largely compensated for by B c — > duj sc + gg where the gluon hadronizes mainly 
into light-flavour hadrons. Details will be discussed in a future publication. 

5 Conclusions 

B c decays represent a particularly intriguing lab to study the interplay of strong and 
weak forces in a non-trivial environment. The I/ttlq expansion derived from QCD 
makes clear predictions on the global pattern: 
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• a short B c lifetime well below 1 psec; 

• a preponderance of charm over beauty decays among the non-leptonic modes; and 

• a reduced semileptonic branching ratio with roughly equal contributions from b — ► 
civ and c — > slv. 

Essential for the analysis is the observation that in a treatment genuinely based on 
QCD there can be no corrections of order I/tuq that have been introduced in purely 
phenomenological models and play a central role there. 

As far as the numerical predictions are concerned, one has to keep an important 
caveat in mind: the weak link in the analysis is the fact that the charm quark mass 
does not provide a parameter that is very large compared to ordinary hadronic scales. 
Thus the l/m c expansion cannot be expected to be quickly convergent. In principle 
it is conceivable that it might actually fail in charm transitions, say through quark- 
hadron duality becoming inoperational there ||. 
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